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HYDROGEN NEUTRAL OUTFLOWING DISKS OF B[E] SUPERGIANTS 

M. Kraus,"'^ M. Borges Fernandes,^ and F. X. de Araiijo^ 



RESUMEN 



o 
o 



> 
(N 

m 
(N 
O 
l> 

o 
o 



ABSTRACT 

The [Oi] line emission of the LMC B[e] supergiant R126 is modeled with an outflowing disk scenario. We find 
that hydrogen in the disk must be ionized by less than 0.1%, meaning that the disk material is predominantly 
neutral. The free-free emission is calculated from the polar wind, and the minimum density contrast between 
disk and polar wind is found to be ^ 10. 
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1. INTRODUCTION & MOTIVATION 
The nature of the B[c] supergiant stars' disks is 
a long-standing problem (see e.g. Kraus & Mirosh- 
nichenko 2006). Recently, Porter (2003) investigated 
the possibility of dust formation in the pole-on seen 
disk of the LMC B[e] supergiant R126 for two dif- 
ferent approaches: an outflowing disk-forming wind, 
and a Keplerian viscous disk. He found that both 
models failed in reproducing the observed dust and 
free- free emission self-consistently. And he suggested 
a way out of this by allowing for substantial alter- 
ation especially of the exponent of the radial density 
distribution from a classical outflow with p ~ 
to a considerably flatter one of r^^'^. Such a flat- 
ter density profile guarantees that at every location 
in the disk a higher density (and hence opacity) is 
maintained, allowing for more efficient dust conden- 
sation and therefore for an enhanced dust emission 
over the disk. Such a modification finally resulted in 
a good fit to the SED of R 126. 

We have taken high- and low-resolution spectra 
of the same object using the FEROS and Boiler & 
Chivens spectrographs with a 2" aperture diameter, 
in order to observe the inner stellar wind material. 
These spectra show astonishingly strong [Oi] emis- 
sion. Due to the nearly equal ionization potentials of 
O and H, this emission must arise in regions where 
the wind material is neutral in hydrogen, i.e. in the 
disk. To test the modified outflowing disk scenario of 
Porter (2003), we took the parameters of his best-fit 
model and calculated the emerging [Oi] line lumi- 
nosities. These turned out to be at least a factor 
of 50 lower than the observed values, meaning that 
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we need a much higher density within the [Oi] line 
forming region than can be provided by Porter's disk 
model. 

The most severe limitation in his model calcu- 
lations was the assumption that, as in the case of 
classical Be stars, the free-free emission arises in the 
high density disk, while contributions from the polar 
wind are negligible. Fitting the near-IR part of the 
SED with free-free emission from the outflowing disk 
therefore determines the disk mass loss rate. Con- 
sequently, this value is an upper limit for the disk 
density and hence opacity, and efficient dust forma- 
tion can severely be hampered. 

2. THE HYDROGEN NEUTRAL DISK 

In our outflowing disk model the density distribu- 
tion follows from the mass continuity equation, i.e. 
p ~ r~^ . The proportionality factor is not limited 
by the free-free emission, but follows purely from the 
fltting of the [Oi] line luminosities. The forbidden 
lines are optically thin. They arise from the five low- 
est levels in Oi that are populated by collisions with 
free electrons. We calculate the level populations 
by solving the statistical equilibrium equations in a 
5-level atom. Collision parameters are taken from 
Mendoza (1983) and atomic parameters from Wiese 
et al. (1966) and from Kafatos & Lynch (1980). The 
oxygen and electron densities are parametrized in 
terms of the hydrogen particle density, i.e. no(r) = 
qon-tiif) and ne{r) = qen}i{r). For qo we use a typ- 
ical LMC abundance value of 1/3 solar, with a solar 
oxygen abundance of go, solar = 6.76 x 10~* (Grevesse 
& Sauval 1998), while follows from the simulta- 
neous fitting of the line ratios and luminosities of 
the [Oi] lines (see Kraus et al. 2007). The electron 
temperature, density parameter (i.e. mass flux over 
terminal velocity), and ionization fraction found for 
the disk- forming wind are given in Table [T] 



1 



2 



KRAUS, BORGES FERNANDES, & DE ARAUJO 



TABLE 1 

DISK AND POLAR WIND PARAMETERS. 



Wind 


Tc 


^m/^oo 


9e 




[K] 


[gcm~3] 




Polar 


10 000 


2.0 X IQ-^^ 


1.0 


Disk-forming 


8 000 


2.2 X 10~" 


4.0 X 10-'' 



To compare our results with those of Porter 
(2003), we need to know the disk mass loss rate. For 
this, the terminal velocity of the disk-forming wind 
needs to be known. The disk of R126 is seen polc- 
on and for the half opening angle wc use Porter's 
value of 10°. Then we can estimate the outflow ve- 
locity from fitting the line profiles of the [Oi] lines 
(see Kraus et al. 2007). A maximum line-of-sight 
component of ^ 2kms"^ might be hidden in the 
line profiles resulting in a disk terminal velocity of 
about ll.Skms""'. The total disk mass loss rate, 
which must be considered as a lower limit, results in 
Mdisk — 2.5 X 10~^ M0yr~^. This is a about a factor 
of 10 higher than the value used by Porter (2003). 

Based on recent observations with the Spitzer 
Space Telescope Infrared Spectrograph, Kastncr et 
al. (2006) derived a total dust mass of the disk 
around R126 of ~ 3 x lO^^Mg. Converting their 
dust mass into a disk mass loss rate we find Mdisk — 
3.4 X 10^"' M0yr~^. This is even 10 times higher than 
our value, confirming the correctness of our derived 
lower limit. 

3. THE POLAR WIND 

With a disk ionization fraction of less than 0.1% 
the disk cannot be the main source of the free-free 
emission. This emission must arise from the polar 
wind. We fitted the SED of R 126 in the optical and 
IR part (Kraus et al. 2007). The best fit polar wind 
parameters are summarized in Tabic [T] Comparing 
the density parameters in the disk and polar wind 
results in a density contrast of ~ 10. This is found 
to be a lower limit. 

4. CONCLUSIONS 

We proposed and tested the scenario of a hydro- 
gen neutral disk for the LMC B[c] supergiant R126 
by modeling the line luminosities and line ratios of 
the [Oi] emission lines resolved in our high-resolution 
optical spectra. The parameters derived for the disk 
and wind of R 126 are the following: 

• We found that the [Oi] 6300 A/5577 A line ratio 
is very sensitive to the ionization fraction in the 
disk (see Kraus et al. 2007). From fitting the 



observed line ratio we can conclude that hydro- 
gen in the disk is ionized by less than 0.1%. This 
confirms that the disk is indeed predominantly 
neutral in hydrogen. 

• The disk mass loss rate of Afdisk ^ 2.5 x 
10~^MQyr~^ found from our fitting is about a 
factor of 10 higher than the value used by Porter 
(2003). This is in good agreement with the pos- 
tulated need for a higher disk density in order 
to fit the [Oi] lines, and also in good agreement 
with the total dust mass of about 3x 10"'^ M0 as 
derived by Kastncr et al. (2006). A disk with 
(much) higher density provides a much better 
environment for efficient dust condensation than 
the disk model used by Porter (2003). 

• The near-IR excess is fitted with free-free and 
free-bound emission from the B-type line-driven 
polar wind. The resulting density contrast be- 
tween equatorial and polar wind is on the order 
of 10. This value is found to be a lower limit. 

To summarize, based on a detailed investigation of 
the emerging [Ol] emission lines we found that the 
disk around the B[e] supergiant R126 is neutral in 
hydrogen right from the stellar surface. Since all 
B[e] supergiants show strong [Oi] line emission, wc 
postulate that the hydrogen neutral ouflowing disk 
scenario might also hold for the other members of 
the B[c] supergiant class. 
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